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The predictability problem in the inverse energy cascade 
of two-dimensional turbulence is addressed by means of di- 
rect numerical simulations. The growth rate as a function of 
the error level is determined by means of a finite size exten- 
sion of the Lyapunov exponent. For error within the inertial 
range, the linear growth of the error energy, predicted by 
dimensional argument, is verified with great accuracy. Our 
numerical findings are in close agreement with the result of 
TFM closure approximation. 

PACS NUMBERS: 47.27Gs, 47.27Eq, 05.45Jn 



Unpredictability is an essential property of turbulent 
flows. Turbulence is characterized by a large number of 
degrees of freedom interacting with a nonlinear dynamics. 
Thus turbulence is chaotic (and hence unpredictable), 
but the standard approach of dynamical system theory is 
not sufficient to characterize predictability in turbulence. 

In fully developed turbulence, the maximum Lyapunov 
exponent diverges with the Reynolds number thus be- 
ing very large for typical turbulent flows. Nevertheless, 
a large value of the Lyapunov exponent do not imply 
automatically shot time predictability. A familiar ex- 
ample is the atmosphere dynamics: Convective motions 
in the atmosphere make the small scale features unpre- 
dictable after one hour or less, but large scale dynamics 
can be predicted for several days, as it is demonstrated 
by weather forecasting. This effect, which can be called 
"strong chaos with weak butterfly effect" arises in sys- 
tems possessing many characteristic scales and times. 
From this point of view, turbulence probably represents 
the example most extensively studied. 

The first attempts to the study of predictability in tur- 
bulence dates back to the pioneering work of Lorenz |l]] 
and to Kraichnan and Leith papers pUS] . On the basis of 
closure approximations, it was possible to obtain quanti- 
tative predictions on the evolution of the error in different 
turbulent situations, both in two and three dimensions. 

A more recent approach to the problem is based on 
dynamical system theory. One of the first results is the 
Ruelle prediction on the scaling of the Lyapunov expo- 
nent with the Reynolds number Q|. Chaotic properties 
have been extensively investigated in simplified models 
of turbulence, called Shell Models, with particular em- 
phasis on the relations with intermittency |||g] . Because 
predictability experiments in fully developed turbulence 
are numerically rather expensive, a similar study on di- 
rect numerical simulations of Navier-Stokes equations is 



still lacking. 

In this letter we address the predictability problem for 
two-dimensional turbulence by means of high resolution 
direct numerical simulations. Turbulence is generated in 
the inverse cascade regime where a robust energy cascade 
is observed Q| . The absence of intermittency corrections 
makes the problem simpler than in the three-dimensional 
case: velocity statistics (energy spectrum, structure func- 
tions) is found to be in close agreement with self-similar 
theory a la Kolmogorov. 




FIG. 1. Stationary energy spectrum E(k) (thick line) and 
error spectrum E&(k,t) at time t = 0.1,0.2,0.4,0.8,1.6. 
kf = 320 is the forcing wavenumber. In the inset we plot 
the compensated spectrum e _2 ' 3 fc 5 ' 3 E(k). 

The model equation is the two-dimensional Navier- 
Stokes equation written for the scalar vorticity co(r,t) = 
—Aip(r, t) with generalized dissipation and linear friction 



d t uj + J(w, i/i) = (-1) p+1 v p A p uj -au)-f 



(1) 



where J represent the Jacobian with the stream function 
-0 from which the velocity is v = (dyip, —d x ip). p is the 
order of the dissipation, p — 1 for ordinary dissipation, 
p > 1 for hyperviscosity. As it is customary in numeri- 
cal simulations, we use hyperviscous dissipation (p = 8) 
in order to extend the inertial range. Although this can 
affect the small scale features of the vorticity field || , in 
our simulations dissipation is not involved in the cascade 
and has simply the role of removing enstrophy at small 
scales. The friction term in ([!]) removes energy at large 
scales: it is necessary in order to avoid Bose-Einstein 
condensation on the gravest mode |10| and to obtain a 
stationary state. Energy is injected into the system by 
a random forcing (5-correlated in time which is active on 
a shell of wavenumbers around kf only. Numerical inte- 
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gration of (|l|) is performed by a standard pseudo-spectral 
code fully dealiased with second-order Adams-Bashforth 
time stepping on a doubly periodic square domain with 
resolution N = 1024. 

Stationary turbulent flow is obtained after a very long 
simulation starting from a zero vorticity initial field. 
At stationarity we observe a wide inertial range with 
a well developed Kolmogorov energy spectrum E{k) = 
Ce 2 / 3 fc~ 5 / 3 (Figure |l|). Structure functions in physical 
space are found in agreement with the self-similar Kol- 
mogorov theory . 

Starting from the stationary configuration, the pre- 
dictability experiment integrates two different realiza- 
tions of the turbulent field and looks at the evolution 
of the difference (or error) field defined for the velocity 
coordinates as 

*u(r,t) = -^(ui(r,t)-u 2 (r,*)) (2) 

From (Q) one defines the error energy and the error 
energy spectrum as ||^] 

E*(t) = J E A {k,t)dk = ^ |<5u(r,i)| 2 d 2 r (3) 

Normalization in (|J) ensures that E A {k,t) — > E(k) for 
uncorrelated fields (i.e. for t — » oo). 

Assuming infinitesimal initial errors, the magnitude 
of the difference field starts to grow exponentially and 
E A (t) ~ E A {Q)exp{2Xt) where A is the maximum Lya- 
punov exponent of the system. The error growth is this 
stage is confined at the faster scales in the inertial range, 
corresponding in our model to the scales close to the forc- 
ing wavenumber kf (see Figure 0). For finite errors, when 
E&(kf, t) becomes comparable with E(kf), the exponen- 
tial growth terminate and an algebraic regime sets in. 
The dimensional prediction proposed by Lorenz |Q as- 
sumes that the time it takes for the error to induce a com- 
plete uncertainty at wavenumber k is proportional to the 
characteristic time at that scale, t ~ T~(k). Within the 
Kolmogorov framework, r(k) ~ e _1 / 3 fc~ 2 / 3 . At larger 
scales the error is still very small in comparison with the 
typical energy. We thus can write 

E A (k',t = r(k)) = [° m (4) 

By inserting (Q) in (^), using the Kolmogorov spectrum 
for E(k) and inverting the dimensional expression for 
r(fc) one ends with the prediction p|,pl| 

E A (t) = Get (5) 

The numerical constant G in (JbJ) can be obtained only 
by repeating the argument more formally within a clo- 
sure framework |^,^,|l2|. In Figure || we plot the time 
evolution of the error energy (E A (t)) obtained from di- 
rect numerical simulations averaged over 20 realizations. 



Both the exponential regime a small time and the linear 
regime (||) are visible. 




5 10 15 20 25 30 35 40 45 50 
t 

FIG. 2. Average error energy (E&(t)) growth. Dashed line 
represents closure prediction (H), dotted line is the saturation 
value E. The initial exponential growth is emphasized by the 
lin-log plot in the inset. 

The dimensional predictability argument given above 
can be rephrased in a language more close to dynamical 
systems by introducing the Finite Size Lyapunov Expo- 
nent analysis. FSLE is a generalization of the Lyapunov 
exponent to finite size errors, which was recently pro- 
posed for the analysis of systems with many characteris- 
tic scales fl3|-|l5||. In a nutshell one computes the "error 
doubling time" T r (S), i.e. the time it takes for an error 
of size S = \Su\ to grow of a factor r (for r — 2 we have 
actually a doubling time). The FSLE is defined in term 
of the average doubling time as 

X{s) = wm^ (6) 

It is easy to show that definition (||) reduces to the stan- 
dard Lyapunov exponent A in the infinitesimal error limit 
@. For finite error, the FSLE measures the effec- 
tive error growth rate at error size 6. Let us remark that 
taking averages at fixed time, as in (^) is not the same of 
averaging at fixed error size, as in (^). This is particularly 
true in the case of intermittent systems, in which strong 
fluctuations of the error in different realizations can hide 
scaling laws like (^|) Jl6[ ] . From a numerical point of view, 
the computation of X(S) is not more expensive than the 
computation of the Lyapunov exponent with a standard 
algorithm. 

The same dimensional argument leading to (|5|), re- 
peated for T r (5), gives the prediction for the FSLE in 
energy cascade inertial range 

X(S) = Aed- 2 (7) 

where the constant A is again not determined by dimen- 
sional arguments. 

The scaling (Q) , which can be shown to be not affected 
by possible intermittency corrections (as in 3D turbu- 
lence [jl3| ), is valid within the inertial range u(kf) < S < 
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U where u(kf) represents the typical velocity fluctuation 
at forcing wavenumber and U ~ y/2E is the large scale 
velocity. 



TURBO", is gratefully acknowledged. Numerical simu- 
lations were partially performed at CINECA within the 
project "Fully developed two-dimensional turbulence". 
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FIG. 3. Finite size Lyapunov exponent X(S) as a function 
of velocity uncertainty S. The asymptotic constant value for 
5 — > is the maximum Lyapunov exponent of the turbulent 
flow. Dashed line represent prediction {]?]). In the inset we 
show in the compensated plot X(8)S 2 /e. The line represent 
the fit to the constant A ~ 3.9. 

Figure || shows the FSLE computed from our simula- 
tions. For small errors S < u(kf) (corresponding to an 
error spectrum E/\{kf,t) « E(kf)) we observe the con- 
vergence of A (S) to the leading Lyapunov exponent. Its 
value is essentially the inverse of the smallest character- 
istic time in the system and represents the growth rate 
of the most unstable features. At larger 6 > 10~ 2 we 
clearly see the transition to the inertial range scaling (0). 
At further large S ~ U ~ 0.1, X(S) falls down to zero in 
correspondence of error saturation. 

In order to emphasize scaling (0), in Figure || we also 
show the compensation of X(S) with e5~ 2 . Prediction (Q) 
is verified with very high accuracy which allows to deter- 
mine the value of A ~ 3.9 ±0.1. The constant A relates 
the energy flux in the cascade to the rate of error growth. 
In absence of intermittency and with r ~ 1 it is possible 
to relate (H) to (0). In the present case (r ~ 1.12) one 
obtains G ~ 4.1. It is interesting to observe that the nu- 
merical result is very close to the old prediction obtained 
by a Test Field Model closure || which gives G = 4.19. 
At least from the point of view of predictability, two- 
dimensional turbulence seems to be very well captured 
by low-order closure scheme. As a consequence we can 
exclude, on the basis of our numerical findings, the ex- 
istence of intermittency effects in the inverse cascade of 
error. This is a result which is probably of more general 
interest than the specific problem discussed in this letter. 
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